Building all Time Evolutions with Rotationally Invariant 

Hamiltonians 



I. Marviaifland R.B. ManiJl] 

Department of Physics, University of Waterloo, 
Waterloo, Ontario, N2L 3G1, Canada 
Institute for Quantum Computing, University of Waterloo, 
Waterloo, Ontario, N2L 3C1, Canada 



Abstract 

All elementary Hamiltonians in nature are expected to be invariant under rotation. Despite this 
restriction, we usually assume that any arbitrary measurement or unitary time evolution can be imple- 
mented on a physical system, an assumption whose validity is not obvious. We introduce two different 
schemes by which any arbitrary time evolution and measurement can be implemented with desired 
accuracy by using rotationally invariant Ifamiltonians that act on the given system and two ancillary 
systems serving as reference frames. These frames specify the z and x directions and are independent 
of the desired time evolution. We also investigate the effects of quantum fluctuations that inevitably 
arise due to usage of a finite system as a reference frame and estimate how fast these fluctuations tend 
to zero when the size of the reference frame tends to infinity. Moreover we prove that for a general 
symmetry any symmetric quantum operations can be implemented just by using symmetric interactions 
and ancillas in the symmetric states. 
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1 Introduction 



In describing any physical system a reference frame is indispensable. The state of a system is not an abstract 
concept independent of a reference frame; the only meaningful and measurable observables are correlations 
between the main physical system and a reference frame. A reference frame itself is always represented with 
a physical system. In practice we do not consider reference frames as quantum systems, and we treat them 
classically. However there are theoretical considerations that force us to treat reference frames quantum 
mechanically. This can be regarded as the counterpart of the famous idea that information is physical: since 
information is stored in the state of some physical object that is subject to the laws of physics, at the most 
fundamental level the abilities and limitations of these physical objects for information processing depends 
on these physical laws. In the same way, a reference frame is always carried by a physical system. Indeed, 
it has been argued that the various conceptual conundrums in quantum mechanics - such as controversies 
about existence of superposition between number states in quantum optics, superconductivity and Bose- 
Einstein condensation, or the problem of quantification of entanglement in systems of bosons and fermions 
[1] - are rooted in ignoring the role of reference frames. 

Amongst the variety of reasons for considering reference frames as quantum mechanical objects, we 
are interested in the following. First, all elementary interactions in nature are expected to have specific 
symmetries. For example, they do not have a preferred direction and so are rotationally invariant. Given 
this situation, a question arises as to whether or not it is possible to implement an arbitrary Hamiltonian 
and measurements that might not have that symmetry. The problem of restrictions on quantum operations 
imposed by assuming a given symmetry was first emphasized by Wigner [2]. Assuming all Hamiltonians 
commute with some symmetry operator, he showed that it is impossible to measure an observable that does 
not commute with that symmetry operator. Later Araki and Yanase [3] proposed a scheme to measure 
an operator that does not commute with that symmetry operator to any desired accuracy. To this end 
they used a quantum system acting as a reference frame beside the main system. The accuracy of that 
measurement increases with the size of the Hilbert space of that reference frame (though not necessarily 
linearly) . 

The second reason for treating reference frames as quantum mechanical objects, which is related to the 
first one, arises from the aspiration to have a relational quantum mechanics [3J[5]- By relational we mean a 
theory in which we do not have an external reference frame that is used to define kinematical observables. 
The key motivator here is general relativity, a background independent theory whose quantum generalization 
is also expected to inherit this property. In the semi-classical limit, where reference frames are assumed to 
have a large Hilbert space, relational quantum mechanics reduces to standard quantum mechanics. However 
at the semi-classical limit we cannot neglect gravitational effects due to these reference frames |H[3J|S1|7]. 
So the problem of finding "relational quantum mechanics" is necessarily part of the project of constructing 
a quantum theory of gravity. 

There are other situations that force us to treat reference frames as quantum mechanical objects. Specif- 
ically, in quantum computation we need to perform measurements on a small quantum register. In this case 
the quantum register must be strongly coupled to the apparatus. This requirement forces usage of small 
and cold apparatus [S]. For example, in some proposed experiments for measuring spin, a small magnet is 
used to measure a single spin. Note that this magnet can be considered as the reference frame that defines 
a specific direction. Due to its small size, the magnet should be treated as a quantum mechanical object 

mm ■ 

Considering a bounded quantum mechanical object as a reference frame has some side effects. First of all. 
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using a bounded quantum reference frame we cannot implement the measurement or time evolution perfectly. 
There are inevitably quantum mechanical fluctuations [T1[51[TT]. Second, performing a measurement or a 
time evolution has always some back-reaction on the reference frame and degrades it (T) EJ [TT] . 

We consider here the problem of implementing any arbitrary unitary transformation on a physical system. 
We show that this can be accomplished using two ancillary systems that act as two reference frames, referred 
to as X-RF and Z-RF, respectively specifying the x and z directions. Note that the state of the reference 
frame is independent of the desired unitary transformation. As the dimension of the Hilbert space of both 
reference frames tends to infinity, the implemented operation becomes the same as the desired one. We 
first construct a primary scheme in which, by using an isotropic unitary transformation and a reference 
frame Z-RF, we can implement all unitary operations on the system that commute with Lz- Wc then use 
this to construct two different schemes, scheme I and scheme II , each of which ensures that our desired 
unitary transformation is implemented to a given accuracy. In scheme I we employ a unitary transformation 
commuting with Lz that acts on both the physical system and an auxiliary frame X-RF whose angular 
momentum is large. Changes in the angular momentum of the system that occur due to the implementation 
of this unitary transformation can be compensated by changes in the angular momentum of X-RF. In the 
scheme II we construct a sequence of unitary operators that alternately commute with Lz and Lx, and 
respectively make use of a sequence of Z-RFs and X-RFs, and show how any desired unitary transformation 
can be implemented via this sequence. 

Any measurement consists of a process in which a measurement apparatus initially uncorrelated with a 
system of interest changes so that some property of the apparatus becomes correlated with some physical 
property of the system. As with all other quantum phenomena, this can be described by a unitary time 
evolution acting on the system and measurement apparatus. Therefore using this scheme we can implement 
any arbitrary measurement just using isotropic interactions; in this sense we can regard this work as a kind 
of generalization of ref. [3 . Our work also generalizes the simple model that Poulin [B] used to show how it 
is possible to construct relational quantum mechanics from ordinary quantum mechanics using some specific 
rules. In that toy model a reference frame specifying the z direction was constructed and used to measure 
the spin of a spin- 1/2 particle. 

Alternatively, this problem is closely related to the issue of restrictions imposed by superselection rules. 
Originally these rules were regarded as some axiomatic limitations added to the theory that restrict the set 
of physically realizable operations [12]. For example the superselection rule for electric charge forbids the 
production of superposed states with different charges. On the other hand, as emphasized more recently [1], 
the lack of a reference frame breaking some symmetry imposes restrictions on implementable operations, 
which can can be regarded as superselection rules. For example, assuming all interactions are rotationally 
invariant, without any reference frame that breaks this symmetry, we cannot prepare states that are not 
rotationally invariant; we also cannot perform time evolutions and measurements that are not rotationally 
invariant. In this way the latter superselection rules might seem more fundamental. However as first argued 
by Aharanov and Susskind |13l I14j , emergence of those fundamental superselection rules in non-relativistic 
quantum theory results from making the (unphysical) assumption that there exist absolute operations 
without any dependence on some quantum reference frame. In ref. |15j Kitaev et al. showed that all 
superselection rules described by compact groups result from the lack of an appropriate reference frame. In 
their terminology, an invariant world, which is subject to a superselection rule described by a compact Lie 
group by virtue of an unbounded quantum reference frame (which itself is subject to those superselection 
rules) can simulate the physics of an unrestricted world that is not subject to those superselection rules. To 
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show this, they prove that for all states, time evolutions, and measurements on a particular system in the 
unrestricted world there exist states, time evolutions and measurements associated with the total system, 
including that particular system and quantum reference frame, that are invariant under the group describing 
those superselection rules. Moreover they produce the same observable effects. Then using this fact they 
show that superselection rules do not enhance the information-theoretic security of quantum cryptographic 
protocols. 

However, our work is distinct from ref. jl5j in that we wish to perform an arbitrary time evolution or 
measurement on some given unknown arbitrary state of a physical system; we do not want to simulate that 
physics, but we want to exactly perform unrestricted operations using a quantum reference frame which 
breaks that symmetry. We therefore assume that the initial states of the system and quantum reference 
frame are uncorrelated. However in the scheme proposed by Kitaev et al [TS], to a state of the system in 
the unrestricted world they associate an initial state in the invariant world that is not necessarily the tensor 
product of that initial state of the system and a state of the reference frame. Indeed, the only case in which 
this initial state is a tensor product is just the trivial case where the initial state of the system is invariant 
under that group. Hence given some unknown arbitrary state of the system, using their scheme we cannot 
perform all quantum operations. 

The preceding discussion was concerned with cases in which implementation of non-symmetric time 
evolution requires non-symmetric resources, i.e. a reference frame that breaks the symmetry. However it is 
usually implicitly assumed that to implement a symmetric quantum operation one needs no non-symmetric 
resources. In the other words, the assumption is that any symmetric quantum operation can be realized by 
a symmetric unitary evolution acting on the system and an ancillary system in a symmetric state. Before 
working more on the non-symmetric cases we check the validity of this assumption in the next section. For 
a general compact symmetry group, and any given symmetric quantum operation we find an explicit form 
of a symmetric unitary time evolution acting on the system and a symmetric ancilla that realize the given 
quantum operation. |17j 

In this paper we denote by i?-inv a rotationally invariant quantity, and Z-inv,X-inv for unitary operators 
that are invariant under rotation around the Z and X axes respectively. 

The outline of our paper is as follows. In section [2] we show that only symmetric resources are required 
to carry out any given symmetric quantum operation. Then in section [3] we introduce a distance measure 
for quantum operations and investigate some of its useful properties. In section [4] we discuss how to classify 
unitary operators with rotational symmetry, both in general and about one axis. We discuss in section |5] 
how to implement Z-inv unitary operators using i?-inv (ie fully rotationally invariant) operators. In section 
[6] we show how to use Z-inv unitaries to implement all unitary operations, and then in the next two sections 
we respectively describe each of schemes I and II. We close in section |9] by comparing the two schemes from 
the point of view of the resources they require and with some suggestions for further research. 

2 When are non-symmetric resources necessary? 

Suppose a unitary time evolution with a specific symmetry acts on the main system and an ancillary system. 
We assume the initial state of ancillary system is also invariant under the symmetry. It is easy to see that 
the total effect of this time evolution on the system is described by a quantum operation that is invariant 
under the symmetry. Clearly if we make use only of symmetric resources, i.e. symmetric time evolutions 
and symmetric initial states, we obtain a symmetric quantum operation. 
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Consider next the inverse problem: is it possible to implement any given symmetric quantum operation 
using only symmetric resources, i.e. with a symmetric unitary time evolution and with an ancillary system 
which initially is in a symmetric state. As we might guess intuitively the answer is yes as we shall now 
demonstrate. 

Consider a quantum operation e that is invariant under some group described by G. We call it a G-inv 
quantum operation. Our aim is to show that for any group G and for any given G-inv quantum operation 
there exists a G-inv unitary time evolution ^5^ acting on the system and ancilla such that 

£(p)=tw(^(p®|0)(0|)^t) (2.1) 

where |0), the initial state of the ancilla, is also G-inv. We prove this fact and give an explicit form of such 
a unitary time evolution |17j . 

We begin by noting that it has been shown |18j that a G-inv quantum operation always admits a Kraus 
decomposition with Kraus operators Kj„ia, where j denotes an irrep, m a basis for the irrep, and a a 
multiplicity index, satisfying 

ng)K,rao.T\g) = ^ V.g € G (2.2) 

where T(g) is the unitary operator in the Hilbert space of the system representing the effect of g € G on the 
state of system and u^-'^ is an irreducible unitary representation of G. We define \j, n, a) a basis in which 
u^^\g) has the following form. 

T{g)\],n,a) = J^')(ff)|j, n, a) ^Y,u^)^{g)\j,n\a) (2.3) 

n' 

On the other hand, if u'-'' is a representation of G, then complex conjugate of m^^-* in any specific basis is 
also a representation of G which might be equivalent to u'^^y We denote the representation obtained by 
complex conjugate of u^^\g) in the basis |j, ri, a) by u^^\ There exists a basis \j,n,a) such that 

T{g)\j;^) = u^^\g)\j;^) = Y.^^^ni9)\i:<^) (2.4) 



To purify e we assume we have an ancillary system initially in the G-inv state |0) as in eq. (2.11 . We 
want to build a G-inv unitary ,5^ acting on the system and an ancillary system such that the effect of this 
unitary time evolution on the system is described by e. Assume jV') is an arbitrary state of system. We 
define Pq to be the subspace spanned by all states of the form of |V') |0) . Also we define Pi to be the subspace 
spanned by all states S'|V')|0) where S is 

sm^) = Kj^^mj;^) (2.5) 



where |j, m, a) are states of ancillary system for which Eq.(2.4l holds. So by definition S* is a map from Pq 
to Pi .We assume |0) is chosen such that it is orthogonal to all states |j, m, a) appearing in Eq.(2.5l. This 
is possible because the ancilla can always be taken to have any number of singlet representations, one of 
which can be taken to be |0). So obviously Pq and Pi are orthogonal to each other. 

Using the normalization condition m a ^Ima^jma = / it is straightforward to see that this map 



preserves inner product. Moreover we can show that S commutes with G. Using eqs. (2.2 1 and (2.4 1 we 
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have 

jma 

= EE"™'™(5)^.™'"n5)i^)®E"n'™(3)i^^) (2-6) 

jma m' n' 

and from the unitarity of u^^"^ wc find 

T{g) (g> T{g) ^| V) |0) = ^ K,,,^T{gm ® |j, n, a) (2.7) 

Therefore 

{T{g)<»T{g)) Sm^) = S{T{g) ^ im\Q) = S{T{g) <^^T{g))mQ) (2.8) 

where in the last step we have used this fact that |0) is G-inv. So S commutes with G. Note that 
T{g) ® T{g) |V')|0) is still in Pq and so T{g) ® T{g)S\i:)\0) is still in Pi. 
Now we define S^^ to be a map from Pi to Pq such that 



V if,™a|V')|j,"i,a) = |^)|0) (2.9) 



By definition all states in Pi are of the form X^jma ^jrna\^)\j, "m, a) for some j?/;) and so S* ^ is defined for 



all states in Pi. Since S preserves inner product so does S . Also from Eq.(2.8l we can deduce that S 



commutes with G. We define the subspace P2 to be the subspace of all states in the Hilbert space of system 
and ancillary system except those who live in Pq find Pi. Finally we can define a G-inv unitary in the 
following manner. 

e Po : = S\^) 

|f^)ePi : .S^]^) ^ S-^\^l) (2.10) 
\n) e P2 : = \n) 

So exchanges Pq by Pi and leave all states in P2 unchanged. Obviously 5^ is unitary and commutes with 



G. Also it is straightforward to check that it satisfies Eq.(2.1l and so its total effect on the system would 
be e. This means that for any G-inv quantum operation one can find a G-inv unitary time evolution such 
that effect of this unitary on the system is described by this quantum operation. 

In general if a quantum operation is described by N independent Kraus operators, to implement it with 
a unitary time evolution we need an ancillary system with an A'^ dimensional Hilbert space. In the scheme 
proposed above for implementing a G-inv quantum operation described by N independent Kraus operators, 
we need an A + 1 dimensional ancillary. So we see that restricting to symmetric resources increase the 
minimum of size required ancillary system at most by one. 

With the same kind of argument we can prove that any given G-inv density operator of system can be 
purified such that the total pure state is G-inv. By definition a G-inv density operator should commute 
with all members of this group. 

Any unitary representation of a group G allows a decomposition of the Hilbert space into charge sectors 
M'j where each charge sector carries an inequivalent representation Tj of G. 



=^-0^- (2.11) 
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Each of these sectors has n{j) copies of the representation j. On each of these sectors the effect of the 
representation of 5 € G, T{g), can be factorized to the irreducible representation associated with j, Tj{g), 
times an identity that acts on the multiplicity subsystem /„(j). Hence it can be written in the form 
Tj{g) (g) /„(j). Each sector can be therefore be decomposed into virtual subsystems [19] 

Mj^^j<g).Aj (2.12) 

The effect of T(g) on .^^j is Tj(g) and on yKj is trivial . So finally the representation T(g) can be written 
in the form 

n5)-0T,(g)®/„o) (2.13) 

j 

We call ^j's gauge spaces and ^'s multiplicity spaces. In the language of quantum information and 
are called decoherence-full subsystems and noiseless . 

Using Schur's lemmas we can show that a G-inv density operator can always be written as 

PG-.nv^Y.P^J^^^P^'^ (2.14) 

where j specifies different inequivalent irreps, is the identity operator on the gauge subsystem of each 
irrep, pj is a density operator acting on the multiplicity subsystem of each irrep, and {pj} is a probability 
distribution. Suppose l^'--'''), which is a purification of p^^\ is expanded as follows: 

1^0)) 1^)1/3) . 

We can deduce that the following state is a G -inv purificcition of pc—inv 

where |j, n) and |j, n) are respectivel vectors in the gauge subsystem of the main system and the ancillary 



one. With the same kind of argument we used in Eq.(2.6) we can see that this state is G-inv. As a simple 



example, we can check that for the case of SU{2) as the symmetry group and the invariant density operator 
for spin half, which is 1/2, the invariant purification state given by this method would be the singlet state. 

Suppose Alice's world is interacting and probably entangled with Bob's world such that all of states 
and time evolution in her world have a specific symmetry. According to the above discussion she can never 
ascertain whether the whole world, including her system and Bob's, has this same symmetry. On the other 
hand, if there existed quantum operations or mixed states that were not purifiable by symmetric resource, 
by observing those states or time evolutions she could ascertain that the whole world is not subject to that 
symmetry. 



3 Distance between two quantum operations 

To estimate the accuracy of implementing unitary time evolutions we need a measure to compare the 
implemented time evolution with the desired one. In this section we introduce a distance for comparing two 
quantum operations and investigate some of its main properties that we shall utilize in this paper. 

Consider two quantum operations ei and £2 that map density matrices to density matrices. By definition 
these are trace-preserving completely positive maps. Consider performing all Von- Neumann measurements 
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on ei{p) and S2{p), and let d{ei,e2) be the maximum difference in probability between the same specified 
outcomes. This comparison can be repeated for different density matrices. We therefore define a measure 



d(£i,e2) =max|ir(^'[ei(p)-£2(p)])| (3.1) 

where the maximization is taken over all density operators p and different outcomes of all measurements that 
can be described by projectors P. This provides a natural measure of the similarity of two superoperators. 
Clearly the distance between two quantum operations varies between zero and one. 

We pause to discuss some of the properties of this distance operator. Consider a unitary time evolution 
Un acting on an N dimensional Hilbert space where iV is a large number. Suppose we partition this N 
dimensional Hilbert space to a pair of subspaces with N — 2 and 2 dimension. Assume the unitary time 
evolution is in the form of In-2®U2 where is the identity operator which acts on the N ~2 dimensional 
subspace. Also U2 acts unitarily on the 2 dimensional subspace. So the unitary time evolution described 
by Un acts the same as the identity operator through a large part of the N dimensional Hilbert space. 
However, the distance of the quantum operation described by ?7Af(-)^iv ^^'^ identity quantum operation 
can be anything between and 1, dependent on the choice of element from U2- Note that an experimentalist 
can easily distinguish these two quantum operations by preparing the initial state of the system to have a 
non- vanishing density matrix only in the 2 dimensional subspace. So two quantum operations with large 



distance (as measured by eq. (3.1)) might act almost the same over a large portion of Hilbert space. Note, 
however, that if the distance of two quantum operations tends to zero, this means that they act almost the 
same in all parts of Hilbert space. This property justifies the usage of this specific distance for our purpose. 
If the distance of a simulated quantum operation with the desired one tends to zero then, no matter what 
is the size of Hilbert space, the simulated quantum operation acts exactly the same as the desired one in all 
parts of Hilbert space and so they are indistinguishable. 

We may also express this definition by making use of the trace distance D(pi,p2) between two density 
operators [TB] 

d{ei,e2) ^ma,xD{ei{p),e2{p)) 
p 

One of the important properties of trace distance is D{e{pi),e{p2)) < D{pi,p2) \l6i- Using this property 
and the triangle inequality we can deduce 

d{e^{e2{.)),e[{s'2{.))) < d{e,{.),e[{.)) + d{62{.),e'2{.)) (3.2) 



Suppose po is the density operator and Po the projector for which this maximum in Eq.(3.1) occurs. Let 
{!«)} be a basis in which ei{po) — £2(^0) is diagonal. Dividing {\i)} into two distinct eigenspaces {|«)^} and 
of positive and negative eigenvalues, we denote the projector to these subspaces by P+ and . The 
trace of ei(po) ~ £2(^0) is zero, so obviously 

MP-[e^{po) - £2(po)])| = MPne^{p„) - £2(^0)])! - ^ J2 " £2(po))N)| (3.3) 



Consequently the projector Pq, for which Eq.(3.1 1 is maximum, can be chosen to be either of P+ or P 
Moreover 



d(ei,£2) = ^ max V|(z|(£i(p)-e2(p))N)| (3.4) 
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where the maximization is taken over all density operators and all orthogonal bases. Using the preceding 
equation we can readily show the convex property 

d{e,^p^e,) <^Ptd{e,e,) (3.5) 

i i 

Suppose e^'^ is a quantum operation acting on systems A and B. Assume initially the state of the 
total system is ® where p^ is some fixed state. The effect of e^'^ on system A is then given by 
e^(p^) = trB{e^'^{p^ ® p^))- Using Eq.(3.4l it is straightforward to check that 



d(£^^(.),£^(.))<d(e^^(.),e^'^(.)) (3.6) 

This inequality expresses the intuitive notion that distinguishing between two different quantum operations 

is more difficult when one observers only part of a larger system. 

The following Lemma (proved in the appendix) provides a useful tool for computing d(ei, £2)- 
Lemma I: Suppose p is a density operator and Oi,02 are two arbitrary operators and {\i)} is an 

arbitrary set of an orthogonal basis. Then 



J2\{i\Oip02\i)\ < IIOill X IIO2II 

i 

Here we have used the infinity norm of an operator, namely 

11^11= max I^J'IV')! (3.7) 

m 

where the maximization is taken over all normalized vectors \tp). Note that ||^i^2|| < ll^ill ll^2||- Using 



Eq.(3.4l and lemma I we can easily see that 



d(c/(.)L/t, y(.)yt) < lit/ _ ^11 +^-\\u- v\\' (3.8) 

So for small ||t/-'t^|| we have d{U{.)W ,V < \\U - V\\. 

4 Classification of unitaries with rotational symmetries 



As we saw in the section 2 , the general form of a unitary representation of a group is given by Eq. (2.131. 
For the rotation group angular momentum is the related index that specifies different representations. So 
the effect of an arbitrary rotation i? on a physical system can be represented by a unitary matrix T(i?) that 
(reducibly) decomposes as 

'max 

T(i?) = TKi?) ® /„(0 (4.1) 

where Ti{R) is the irreducible representation (irrep) of a rotation with angular momentum I. Here n(l) is 
the multiplicity of this angular momentum and /„(;■) is the n dimensional identity. We say that a unitary 
transformation is rotationally invariant iff 

[U,T{R)]^0 (4.2) 

for all rotations R. We denote rotationally invariant unitary operators by U^-inv Noting that Ti{R) is an 
irrep of the rotation group, we deduce using Schur's lemmas that all possible rotationally invariant unitaries 
have the form 



-inv — 
I 



/; C/(') (4.3) 
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Figure 1: Adding a largo vector to a small one: The length of the sums. S and S' of the large vector A and 
the small vectors B and B' is almost independent of the components of B and B' which is perpendicular 
to A. Lemma II demonstrates the quantum version of this simple property. 

where [/*•'•' is an arbitrary unitary operator acting on the n(i)-dimensional Hilbert space of multiplicity 
subsystems and /; is the identity operator on the gauge subsystems. 

Now consider a state \j,m,X) with angular momentum j, magnetic quantum number m, and A some 
other possible quantum number. The effect of J7^_inv on this state is 

«7^-inv|i, m, A) = ^ Ul^,l\j, m, A') (4.4) 

A' 

where Uy\ ^^^^ unitarily on the subspace of A-multiplets with the same j and m. 

Consider next all unitary operators V that have one axis of symmetry. Without loss of generality we 
can take this to be the 2;-axis, ie. we consider all unitaries that commute with L^. They can be decomposed 
via 

V'ir-i„v = VW= E yiy|M,A)(M,A'| (4.5) 

M M,A,A' 

where {|M, A)} for different A is an orthogonal basis for the subspace where Lz = M. The operator F(^) 
acts unitarily on this subspace. 

5 Implementing Z-inv unitaries using R -inv unitaries 

In this section, for any given unitary time evolution with one axis of symmetry, a Z-inv unitary, we construct 
a rotationally invariant unitary time evolution such that, when this R-inv unitary acts on the combined 

physical system with Z-RF, the total effect on the physical system is (xiuivalcnt to that of the given Z-inv 
unitary. As noted in the introduction we refer to this procedure as implementation of a Z-inv unitary. 

The main idea behind this construction is based on the simple property that when two vectors of unequal 
norm are added together, the length of the resultant vector is almost independent of the components of the 
vector of smaller norm that are orthogonal to the larger vector, provided the former has sufficiently small 
norm (see fig. 1). Lemma II demonstrates the quantum version of this simple property. 
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Lemma II: Consider the following state 



= \li,mi)\l2,m2 ^ h - k) (5.1) 

where li,l2 are angular momenta and mi,m2 are eigenvalues of Lz- At the limit of I2 S> li,k^ this state is 

almost the same as the state with total angular momentum equal to mi + 12 and total equal to mi + l2 — k 
i.e. 

10) = \{j = mi + l2,m = mi + I2 - k) : {l^h)) (5.2) 

or more precisely 

> 1 - + (5.3) 

212 

In terms of Clebsch-Gordon coefficients this means 

2 (1? 4-1, - m? 



_ ^ _ {li + li-mi){2k+l) -mi 
2I2 



This lemma is proven in the appendix. A specific case of these relations is reported in [6] where it is 
verified numerically. 
For fc = we have 

l(0l^>P>l-^^i^^^>l-C^ (5.4) 

where — '^^2^;^^^*^ is the minimum of 1 — '^^'^ 2™^ "'^ ' "^tiich occurs for mi = —1/2. 

Now consider an ancillary system called Z-RF with angular momentum Ij^z that is large in comparison 
with the maximum angular momentum of system li. We assume initially Z-RF is in the state |Z — RF) = 
IIrZiHIrz = Irz) S'lid so has maximal angular momentum in the z direction. Following the proof in lemma 
II, to a very good approximation the total angular momentum of the combined system with Z-RF depends 
only on mi (the z component of angular momentum of the system) 

\h,mi,d)\Z - RF) « \ j ^ mi + ljiz,m^ mi + Irz,^) (5.5) 

and so is almost independent of li. Here S labels possible degeneracies of the state of the physical system 
and A labels all distinct states of the total system with the same j — mi + Ifjz and m — mi + Irz- More 
precisely A indicates from which (Zi, Wij^, Z/j^, (5) the state \j = mi + iRZyTn = frii + Irz,^) is formed. 
Since m^z and Irz are fixed, different A's stand for different (^i,(5)'s. 

The upshot is that j and m of the combined system are independent of {li,S), and depend only on 
mi. Hence we can implement any arbitrary unitary on the subspace of states with the same mi by using 
an appropriate U^-inv In other words we can implement all unitary time evolutions on the system that 
commute with Lz acting on the system. In the language of quantum information, we implement an arbitrary 
Z-inv unitary on the system by performing an i?-inv unitary on the noiseless subsystem of the total system, 
which consists of the main system and reference frame. 

Suppose we are going to implement a unitary time evolution operator V on the system that commutes 



with Lz. V can be written in the form of Eq.(4.5 1 



Now consider the following decomposition of total Hilbert space of system and Z-RF induced by rotational 
group 

J^sys «) J^fz-RF ^0 -^3® (5-6) 
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We define the unitary acting on the system and Z-RF to be 



(5.7) 



where Ij acts on the gauge subsystem and ^^j-^^, multiplicity subsystem jVj. Obviously 

"V has the same form as (4.3 1 and so is isotropic. Note that here we assume this unitary governs the total 
system (ie. the physical system combined with Z-RF), and so the context of the above formula is for the 
total system. We assume this unitary acts on the initial state of the total system, which is a tensor product 
of the initial state of the system and \Z — RF). 

The effect of this time evolution on the system is described by a superoperator called e'^^\ One possible 
representation of this time superoperator can be specified by the following set of Kraus operators 



(5.8) 



(5.9) 



Kn = {lHz,n\r\Z-RF) - {lRz,n\r\lRz,lR.z) 
where \lRz,n) is the eigenvector of Lz with eigenvalue n. For if;^^ have 

Ki,, = (iRzjRzinRzjRz) = J2 yx'y""\iRZ,lRz\j,J,X){j,J,X'\lRZ,lRz) 

+ E Vx^''''"\lRzJRz\.hM,X}{j,M,X'\lRzjRz) 
We can rewrite the second term as 

V = {Irz, lRz\{y - - ^)\Irz, Irz) 

where ,^ is the identity operator in the Hilbert space of the combined system and RF, and is the projector 
on the space of all vectors with M = j. So the operator — ^ is the projector to the subspace of states 
with M < j. For an arbitrary normalized vector |0) in the Hilbert space of the physical system we know 
from Eq.(5.4) that \{y - rz A rz)\ < C where = Since r is unitary we deduce that 

IIV^'IKC^. 



To compute the first term in Eq.(5.9) we note that 



(j, j, A|(|/i,m,(5) (g) \IrZ,Irz)) = Sj^m+lRz^^,{h,S) X Cm,A 



(5.10) 



where ^m.x is some real number. From Eq.(5.4l we deduce that ^H, n > 1 — and so , is close to one. 



So using Eq.(5.10l we have 



{lRZ,lRz\j,j, A) = ^j-iaz:\ I'm^j- A) 



(5.11) 



Note that A specifies li and 6. Now using this equality we see that the first term in Eq.(5.9 1 can be rewritten 
as 

Y.^x^''"^^ilRZ,lRz\j,J,X){j,jA'\lRzjRz)^ E Vx7kmMm,x'\m,X){m,X'\ (5.12) 

j.X,X' m,A,A' 



To compute this, first we define 



X = {Irz, lRz\y - ^\Irz, Irz) = ^ - E ^) ^1 

m,A 



(5.13) 
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To get the last equality we have used Eq.(5.11 1. Note that > 1 — and so ||X|| < . Since f„j_A is 
close to one we can easily see that ^m,x ~ 1/2(1 + Cm a)- ^'^'^ Eq. (5.121 equals 

(/ - _ |) y - + xv) (5.14) 

So finally Ki^^^ = V + V where V = V — ^{XV + VX). By the triangle inequality we can see that 
||F|| < 2C2. Finally we obtain 

s^'\p)^VpV^ + ei'l^ip) (5.15) 

where 

^nLeip)^iVpV^ + VpV^)+ J2 ^npKj (5.16) 

n<lRZ 

Note that e^^^^^{p) is not necessarily a positive operator. The first two terms in sli^i^^ip) commute with Lz 
and are responsible for the noise in each block of the same m. The effect of the last term is that of mixing 
states with different to. 

While our target was implementing a unitary time evolution described by V, we have instead imple- 
mented a time evolution that is described by e^^^. To estimate the quality of this implementation we 
compute the distance of these two time evolutions, d{V{.)V^ , e'^^^.)): 

d{Vi.)V^,e('H.))^]-me,xJ2mei%,ip)\^^ J2 K,,pKJm (5.17) 

where the maximization is over all bases {\i)} and all density operators p. To calculate this quantity first 
we note that X]n<;fiz ^npKn is a positive operator and so 



= ^-HKirzPKirJ) 

= l-tr{{V + V)p{V+V)^) 

= -tr{VpV^ + VpV^ + VpV^) (5.18) 

For any set of operators {Oi] we have \tr{Oi . . . Oisip)\ < ||Oi|| . • . ||Ojv||. Using this fact and noting that 
||y|| < 2C^ we can show that for small we have 

tr{ Y KnpKj) = |tr(l/pF^ + Fp0 + FpF^)| < (5.19) 

n<lRz 

On the other hand, we have 

Ymvpv^ + vpv^m < J2 mvpv^\^)\ + E \(^\ypv^m (5-20) 

i i i 

Using lemma I we have 

YmVpV^ + VpV^m\<AC^ (5.21) 

i 

and so 

ci(y(.)Ft, £«(.)) = ^max(^ |(z|eil,(p)|*)|) < (5.22) 
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Note that to get the same amount of error for different systems, the dimension of RF should increase 
proportional to the square of the angular momentum of system. 

After using the reference frame its state would be changed due to back reaction effects. Since the total 
Lz is conserved, if there were no error in the implementation of a Z-mv unitary, the state of the reference 
frame would be unchanged. The probability of error is less than 4C^, so with the probability of 1 — 4C^ or 
more the reference frame stays in its initial state. To get a measure of how the reference frame degrades, 
consider using it to implement a unitary operation on a different system each time |20j. After using the 
reference frame n times, as long as nC^ ^ 1 the probability of being in its initial state would be more 
than 1 — 4nC^ and so the error would be less than (1 — 4nC^)AC^ . We can see that the state of the 
reference frame after the first use would be a mixture of states from the set mf;z — Irz)j \Ir.z, ^tt-rz = 

Irz - 1), \Irz, muz = Irz - 2^i)}. 

For a reference frame with fixed Irz, \Irz, Irz) seems to be the best choice. However, taking into account 
lemma II and the property of adding vectors that we have used in this scheme, we can use other states of 
the form m^^^ — Ijiz ~ k) provided Irz ^ k^- It is straightforward to see in this case instead of 4C^ 
the error would be almost 4C^(2fc + 1). 



6 Implementing arbitrary unitaries using Z-inv unitaries 

In this section we show how to implement all arbitrary unitary time evolutions on a physical system, whether 
or not they commute with L^, by using a Z-inv unitary time evolution acting over that physical system 
combined with an auxiliary reference frame X-RF. We assume the angular momentum of X-RF, Irx, is 
large compared to the maximum angular momentum of the physical system under consideration, Igys- The 
idea is that changes in the angular momentum of the system caused by unitary transformations on the 
system are compensated for by making changes in X-RF. We choose the initial state of X-RF to have a 
large uncertainty in Lz, such that increasing or decreasing its Lz leaves it nearly unchanged. 
We therefore take the initial state of X-RF to be an equal superposition of the form 

\Rx,0)^^^^={\m=.-N) + ...+ \m^N)) (6.1) 

where N — Irx — '^hys, so that the sum is the subset of magnetic quantum numbers within this range. This 
will afford us freedom to compensate for changes in L^ of the system by modifying the state \Rx, 0), which 
we shall denote by \Rx)- Note that the expectation values of Lz and Ly are zero for this state, whereas 
the expectation value of Lx is nonzero, which means that this state is pointing in the X direction. Unlike 
\Z — RF), which is an eigenvector of Lz, \Rx) is not an eigenvector of Lx and has completely different 
character. We shall also deploy the following notation 

\Rx,n) = ^=^==(|m ^-N + n) + ...+ \m = N + n)) (6.2) 

which we refer to as a shift of \Rx) by n. 

Initially X-RF is in the state \Rx) which has a large amount of uncertainty in Lz- To implement 
unitaries on the system that do not commute with Lz, we compensate for the change in Lz by shifting 
the state of X-RF such that the total Lz remains constant. Under this change the final state of X-RF is 
dependent on the state of the system, thereby entangling them. However if N is sufficiently large, which 
means the initial state has a large uncertainty in Lz, these shifts change \X — RF) by only a small amount. 
For very large N the state of system remains uncntanglcd with X-RF. 
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Our goal is to implement any arbitrary unitary transformation U where 

U\m,X) = J2 Uim',X'Um,x)\m',X') (6.3) 

m',A' 

and where |m, A) is a state of the system. To implement U we construct the following Z-inv unitary 

^ = ^^(M) (6.4) 

M 

in which each ^(^) acts in the following manner 

'^(^V,A)|M-m) = ^ C/(„',A'),(m,A)|m',A')|M-m') if \M\ < N + l,y, 

m',A' 

= \m,X)\M-m) if \M\>N + l^y^ (6.5) 

where |m, A) is in the Hilbert space of the main system and \M — m) is in the Hilbert space of X-RF. Using 
the unitarity of U it is straightforward to check that is also a unitary operator that commutes with Lz- 
Now the effect of U on the initial state |m, X)\Rx) is 

^ N+m 

^\m,X)\Rx) = ^^""^HlM-m) 

+ M=-N+m 



^ N+m 

= /oAT+T ^ T.^irn',X'Um,x)\m',X')\M-m') 

M=-N+m m',X' 
^ N+m 

= /oW^ E Uim',X'Um,x)\m',X') Yl \M-m') 

m',X' M=-N+m 

= X] '^('»'.^'),(™.A)l"i', A')|iix,m- m') (6.6) 

For an arbitrary state of the system |9) = J2m X ^m,x\'m, A) we have 

^\e)\Rx) = Y,^rnAYl U^m',x'Um,x)\m',X')\Rx,m-m')) (6.7) 

m,A m',A' 

So after time evolution the state of system becomes entangled with the state of the reference frame. This 
entangled state includes a superposition of vectors that are the tensor product of some state of the system 
and IRxjTn — m'), which is the initial state of the reference frame shifted by m, — m'. Since —Igys < m < lays 
the maximum absolute value of these shifts is 2lays, yielding the extremal shifted states \Rx, ^'^hys) and 
\Rx,+'^lsys)- We are interested in the largest common part between all of these shifted states. Therefore 
we define the unnormalized vector |r) such that the overlap {r\Rx,m — m') is maximal and independent 
of {m, m'} for all allowed values of {m, m'} (see fig. 2). We can easily see that 

N-2ls 



i=-N+2ls 

By writing 



\Rx,m- m') = \T) + {\Rx,m - m') - \T)) (6.9) 
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Figure 2: a) The state \Rx), b) the state |i?x,+2^sys) which is the maximum shifted to the right of \Rx), 
c) the state \Rx^ —'^hys) which is the maximum shifted to the left of \Rx) , d) the unnormalized vector |r) 
which is the common part of all of these states. 



we decompose \Rx,m — m') as superposition of |r) and another vector orthogonal to |r), a result that 
follows from noting that 

{T\Rx,m-m) = {T\T) = 2N + \ 

Using this decomposition we obtain 

'^mRx) = U\Q) ®\V) + Y,e.^.^ Yl U,^m',x'),(m,X)\m\X')i\Rx,m - m') ~ \T)) (6.10) 

m.X m'.A' 

This describes the total state of the system and X-RF after this time evolution. To find the effect of this 
time evolution on the system we trace over the Hilbert space of X-RF. This yields a superoperator that 
maps the initial state of system, p to 

e^'\p) = {V\V)UpU^ + (1 - (r|r))eil,(p) = ^Sll_lllA±lupU^ + (6.11) 



(2) 

where £noise is ^ trace-preserving completely positive super-operator. So using Eq.(3.4l we see that the 
error in the outcome probability of any arbitrary measurement is less than d{U {.)U\ e^'^^) < 2n^+i ■ In the 
limit of large N/lgys the error rate, 2lsys/N, is small. 



7 The First Scheme 

In section [5] we saw how an i?-inv time evolution acting on the physical system and Z-RF can be used to 
implement any arbitrary Z-iiw unitary on the physical system. Then in section[6]we saw how a Z-inv unitary 
time evolution acting on the main system and X-RF can be used to implement any arbitrary unitary on the 
physical system. By combining these two schemes we can perform any arbitrary unitary on the system just 
by using i?-inv unitary time evolutions and Z-RF and X-RF (see fig. 3). 

To implement an arbitrary unitary U acting on the system, we need to implement the Z-inv unitary ^ 
(defined in the previous section) on the system coupled to the auxiliary X-RF. To do so we let the i?-inv 
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All Unitaries 



Figure 3: Schematic description of scheme 1. 

unitary 'V — ^jlj 'rz) act on the system coupled to both X-RF and Z-RF. Initially the density 

matrix is 

p''''(S>\Rx){Rx\(8>\Z-RF){Z-RF\ (7.1) 
After the time evolution the state of the system is 

e{p) = trz-RF,x-RF{yp'''' ® \Rx){Rx\ ^ \Z ~ RF){Z - RF\r^) (7.2) 
is designed to implement on the physical system and X-RF. So 

e{p) = trx-RFi'^p''"' ® \Rx){Rx\^^) + trx-RF{e^^Le{p"''" ® \Rx){Rx\)) (7.3) 



From Eq.(6.11 ) we know that for large N 



e{p) = ^^^UpU^ + ^e,!fi,(p) + trx-RFiei'Lip'^' ^ \Rx){Rx\)) 



^^L^^UpU^ + Sr^^Up) (7.4) 

where 

eno^se{p) = ^ei'ie(p) + ^^X-fli. (^i'le (P^^^ ® \Rx){Rx\)) (7.5) 

To estimate the accuracy of this implementation we find an upper bound for d{U{.)U\e{.)) 
d{U{.)U\e{.)) = imax V|(.|[/p[/t-e(p)|z)| 

i 

- trx-RFie'^^Leip'^' ® \Rx){Rx\m (7-6) 
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Figure 4: Schematic description of scheme 2: building all unitaries using a sequence oi Z — inv and X 
unitaries. 



- mv 



(2) 

Using the triangle inequality and noting that £„ojsg is positive and trace-preserving, we find 
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d(C/(.)C/^£(.)) < ^'^'^ 
21 

< 



N 



1 ^ , 

max 



iV Z p.{\i)}'^ 



Z p,{|i)} ^ 



,{p'y'^\Rx){Rx\)\sM 



(7.7) 



where {\s)} is orthogonal basis in the Hilbert space of X-RF. Using Eq.(5.22) we can find an upper bound 
for the last term in the right-hand side of the above equation. Since the angular momentum of X-RF is 
N + 21 sys, the maximum total angular momentu m of th e system and X-RF is li = N + 3lsys. However since 



N ^ Isys we can assume li m N. Thus using Eq.(5.22l and noting that = '■1+^^+^/'^ -^g obtain 



d{U{.)U\e{.))< 



21 



sys 



2N^ 



N 



RZ 



(7.8) 



where l^z 3> N . Obviously by choosing larger Ijiz the error becomes smaller. For any given Ijiz there is a 
specific N that minimizes this error. After minimizing with respect to N we have 



error i{lsys,lRF) 



2P 

'^'■sys 



1/3 



RZ 



(7.9) 



where error j shows the error in scheme I as a function of the maximum angular momentum l^ys of the 
system and the angular momentum lup of reference frame. We shall discuss the implications of this result 
in section [9l 

8 The Second Scheme 



In this section we propose an alternative method for building an arbitrary unitary time evolution by using 
i?-inv time evolutions. As we have seen in section [s] we can implement all Z-inv unitary evolutions by using 
i?-inv unitary time evolutions and a reference frame in the z direction. In the same way we can build all 
unitaries commuting with Lx (X-inv unitaries) by using a reference frame X-RF that is defined in the 
same way as Z-RF, but rotated so that the x direction is now the specified direction. Note that the X-RF 
we use in this scheme differs from the one we used in the first scheme. Now by a sequence of Z-inv and 
X-inv unitaries we can build more unitary time evolutions. In fact, as we will show in the following, from 
a sequence of unitaries alternately commuting with Lx and Lz we can build any arbitrary unitary (see fig. 
4). 
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The main idea is based on the following property using the four Hamiltonians ±Hi and ±H2- By 
sequentially applying these Hamiltonians we can construct any arbitrary Hamiltonian contained in the Lie 
algebra generated by ^"2}. Consider the following example. Apply Hi, followed hy H^—Hi, and —H2, 
each for the same time 5t. Since 

^iHi5t^iH2St^-tHiSt^-~iH25t ^ ^iiHiH2-H2Hi)St^ _|_ Q(-^^3-j j^-j 

for small St, the result is the same as if one had applied i[Hi,H2] for time St^. In general, any given 
Hamiltonian in the Lie algebra generated by Hi and H2, can effectively be constructed by applying a 
sequence of ±Hi and ±H2 [H]. 

Using i?-inv interactions acting on the system and reference frames, we are able to apply all Hamiltonians 
commuting with Lz and also all Hamiltonians commuting with Lx by using Z-RF and X-RF respectively . 
The Lie algebra generated by these generators describes the full set of Hamiltonians that can be constructed. 
The following lemma illustrates how to find this Lie algebra. 

Lemma III Suppose A,B are two Hermitian operators with the property that no eigensubspace of A is 
orthogonal to any eigensubspace of B. Then the union of the set of all unitaries commuting with A and the 
set of all unitaries commuting with B is a universal set i.e. all unitary operators can be constructed from a 
sequence of unitaries in those sets. Moreover the length of required sequence is uniformly bounded. 

This lemma is proven in the appendix. (It has a nice generalization based on graph connectivity to 
an arbitrary number of operators (instead of just A and B) \22\.) In any representation Lz and Lx have 
the property that no eigensubspace of Lx is orthogonal to an eigensubspace of Lz. Consequently unitary 
operators commuting with Lz and unitary operators commuting with Lx form a universal set. 

Each time we use a reference frame it experiences some inevitable backreaction . For instance after 
performing a Z-inv time evolution on the system, the Z-RF is not in its initial state \Z — RF) — \Irz, Irz) 
anymore, but instead is in a mixture of states including states with other magnetic quantum numbers. This 
new state of Z-RF cannot specify the z direction as well as the initial state. So it cannot be corrected to 
get the initial state of Z-RF by just using R-inv resources without using another Z-RF. We might employ 
this used Z-RF to perform the next Z-inv time evolution; however since this used Z-RF cannot specify z 
direction as well as the initial Z-RF we expect more noise in implementing the second Z-inv time evolution. 
Consequently to avoid increasing amounts of noise we need a fresh reference frame after each use (or after 
several uses) when the Z-RF is degraded more than some specific threshold, dependent on the desired 
accuracy. Considering the reference frames as a resource it is important to know the number needed to 
implement a unitary time evolution with some precision. The length of the required sequence is uniformly 
bounded, which means that there exists an upper bound for the number of steps required to build any 
arbitrary unitary. How can we estimate this number? In other words, how many times do we need to use 
those reference frames? 

To answer this question we employ the Solovay-Kitaev theorem. According to this theorem if G is a 
universal set of unitary operators which produce a dense subset of SU{d), and G is closed under inverse 
then yU e SU{d),T> 0, 3C/i, f7„ G G such that \\U - Ui...Un\ \ < t and n = [16 . 

In the present case G consists of unitaries commuting with Lz and unitaries commuting with Lx. 
Suppose Uapx is the approximation of U that is obtained by this method after n steps Uapx = Ui...Un where 
\\U — UapxW < T. According to the Solovay-Kitaev theorem we can assume n « ^ln^(^) where A is a 
constant. Using Eq.(3.8l and assuming r is small we find d{U{.)U\ Uapxi Wlp^) < T. 

Due to the finite size of the reference frame we cannot perform Ui, ...,Un perfectly; each time there is 
an error less than 4C^. So instead of Ui...Un{.)Ul...Ul we implement e(.) = £i(. ..£„(.)). Using the triangle 
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Figure 5: Error in scheme II as a function of x where x cx iRp/lsys 



inequahty we have 

d([/(0;7t,£(0) <d(C/(0[/^C/ap.(OC^aU) + ^i(C/ap.(Of/aU,£(0) <^ + 4nC2 = T + 4AC2ln2(i) (8.2) 



where Eq. (3.2 1 was repeatedly used to bound d{Uapx{-)Ul , e{.)) . 



The angular momentum of the reference frame, which is proportional to and the required number 

of reference frames n both can be regarded as limiting resources. In terms of these resources total error is 
bounded by 

e-V^ + inC^ 

Obviously increasing the angular momentum of the reference frame decreases and with it the error. As 
— > this error approaches zero. However given a fixed C^, what is the minimum total error we can 
achieve under this restriction? By minimizing the total error t + 4A In^(i) with respect to r we can find 
the minimum accessible error under this constraint. Defining r = e^'^ and x = (8AC'^)~^ the minimizing 
condition is 

X ^ ae" =^ a ^ W{x) (8.3) 
where W{x) is the Lambert-W function [23^. Hence the total error is bounded by 



(see fig. 5) and the number of reference frames we need would be AW'^{x). Note that x is equal to some 

2 

sys 



constant times lup/llys- For a; > 1 we can see that ln(x)/2 < W{x) < ln(a;). So at the limit of large x the 
error is bounded by 



err or J I {I sys, hp) = ln^(x)/a; 
and the number of required reference frames grows like lii^(x). 



(8.4) 
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Finally we note that, though non-optimal, we can of course choose any other direction instead of the 
X direction (except Z) for the second reference frame. However a bad choice of independent direction 
necessitates a longer sequence to reach an arbitrary unitary. 



9 Discussion 

Associated with any restriction on resources is a theory describing whether under certain constraints a 
given operation is feasible. Entanglement resource theory is a well-known example. Analogously, a resource 
theory has more recently been developed for reference frames [18]. It is interesting to compare the two 
schemes proposed here from the viewpoint of the resources they require. 

First note that the error in both of these schemes is a function of Igys/lRF- Hence if the maximum 
angular momentum of the system increases by a factor of k then an increase in the angular momentum of 
the reference frame by a factor of fc^ will yield the same error. It is not clear whether there exist schemes 
in which this factor is smaller than . 

Now let us check how the error decreases as the (large) reference-frame angular momentum increases in 
these two schemes. Looking to eqs. ( 7.9|8.4 1 we see that in the first scheme the reciprocal of the error grows 



as ^{-s/Irf) while in the second scheme it grows as (ffilRp / ^"^^{Irf)) ■ Recall that both these schemes are 
based on a simpler scheme to implement Z — inv unitary time evolutions in which the reciprocal error is 
6'{Irf)- So it seems reasonable that the reciprocal of the error in both of these schemes grows slower than 
^(Irf)- a natural open problem is to determine the best scheme for ensuring the error tends to zero as 
rapidly as possible relative to Irf- 

We have shown that the presence of rotational symmetry does not restrict us in performing arbitrary 
measurements or time evolutions of the system. As there is evidently nothing specific in rotational symmetry, 
we expect that the same holds for the existence of other symmetries. For example an extension of these 
results to more complicated Lie groups, such as SU{N), would be interesting. Here the challenge is to 
construct an appropriate generalization of lemma II. Apart from this, there is nothing specific to SU(2) that 
would appear to prohibit this kind of generalization. Work on this issue is in progress. 
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Appendix: Proofs of the Lemmas 
Lemma I: 

Suppose p is a density operator and Oi,02 are two arbitrary operators and {\i)} is an arbitrary set of an 
orthogonal basis. Then 

Y,\{i\Oip02\i)\ < IIOill X IIO2II 
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Proof 

Using the Cauchy-Schwartz inequality we have 



^ Ki|0ip02|i)| = Y.mO,^p^p02\i)\ < lY,{i\OipOl\i)lj2{i\02pOl\i) (9.1) 

i i V * V * 

= ^tr{\0i\^p)tr{\02\^p) < \\Oi\\ X IIO2II (9.2) 
Using this result we can get the more general result that 

J2\{i\Oi...OhpOk+i...ON\i)\ < ||Oi||...||0;v|| (9.3) 

i 

Lemma II 

Consider the following state 

\ip) = \li,mi)\l2,m2 = l2-k) (9.4) 

where li.l2 arc angular monicntums and mi,m2 arc eigenvalues of Lz- At the limit of I2 3> If^k'^ this 
state is almost the same as the state with total angular momentum equal to mi + 12 and total equal to 
mi + I2 — k i.e. 

10) = |(j = mi + I2, m = mi + l2-k): (li, ^2)) (9.5) 



Zt2 



or more precisely 

|2 > ^ _ {ll + h-ml){2k + l)-mi 
In terms of Clebsch-Gordon coefficients this means 

lim frmi+l2-k,m,+l2\^ _-. {if + h - ml){2k + 1) - mi 

Proof 

Define mtot = mi + I2 — k. For i > mtot we denote 

\i) = \{j = i,m = mtot) ■■ {luk)) (9.7) 
Using this notation we can expand \ip) as 

where Itot = mi+l2 and |(/)), which is defined in the lemma is actually the state |(j = Itot, m — mtot) '■ {h, h))- 
Note that izlV') = "T'totiV') and so vector states with total angular momentum less than mtot do not appear 
in this expansion. Using this expansion we have 

J^\tP) = akot{kot + m)+ Yl + (9-9) 

h+h h+h 

= altotiltot + m}+ E hotiltot + ^)m+ E [ii^ + ^)-ltot{ltot + m^\^) (9.10) 
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h+h 



= itotikot + im + J2 + ^) - hotikot + iMii) 

On the other hand, using the relation 



-2 , r2 



2+ 



(9.11) 



(9.12) 



we obtain 



J2|^) ^ [/^(;^ + 1) + i^^^i^ + 1) + 2mi(/2 - fc)]|^) + 



v/(/i + mi + l)(;i - mi) ^(2^2 - fc)(l + k)\i)i) + - mi + l)(;i + mi)^k{2l2 ~k+ 1)1^-2^) (9.13) 
Where 

IV-i^) = |/i,mi + 1) (g) |?2,m2 = ?2 - fc - 1) , ^ \h-mi - 1) (g) 1^2, ma ^ h - k + 1) 
are orthogonal to 1-0). 



Equating this result with Eq.(9.11 1 we obtain 

h+h 



J2 w + 1) - itotikot + i)]m = Ai^p) + + D\^pi) 



(9.14) 



where 



A = Zi(Zi + 1) - mi(m,i + 1) - 2mifc B = {h + mi + l)(/i - mi)(2;2 -k){l + k) 
D= yj{h - mi + l)(/i + mi)k{2l2 - k + 1) 

So we deduce that 

h+h 

l2 



^ [til + 1) - kotikot + ^A^ + B^ + D' 



(9.15) 



(9.16) 



Now we find the lower bound for the left-hand side of Eq.(9.16 1 for a fixed value of lAP = 1 ~ ct^- The 
minimum of this expression occurs when all (3i are zero except that one for which the factor + 1) — 
hotihot + 1)]^ is minimum. For k > this minimum occurs when i = Itot — 1- For fc = 0, i is larger 



than Itot and so this minimum happens for i = Itot + 1- Putting these /3iS in Eq.(9.16l and noting that 
lap = 1 - E lAP we obtain 



h+h 

A^ + B^ + D^= J2 W+'^)-hot{kot + l)]W>a~\a?mt 



or equivalently 



^kotV 



This inequality is always true for all li, I2, k. Now we assume I2 S> If, k and so we obtain 



2 ^ HI + h - m,l){2k + I) ~ mi 

\a\ > i 

' ' - 2^2 



(9.17) 



(9.18) 



(9.19) 
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Lemma III 

Suppose A,B are two Hermitian operators with the property that no eigensubspace of A is orthogonal to 
any eigensubspace of B. Then the union of the set of all unitaries commuting with A and the set of all 
unitaries commuting with B is a, universal set i.e. all unitary operators can be constructed from a sequence 
of unitaries in those sets. Moreover the length of the required sequence is uniformly bounded. 
Proof 

We define to be the linear space spanned by all Hamiltonians that can be constructed by a sequence 
of applying Hamiltonians commuting with A and Hamiltonians commuting with B. Suppose {|ai, (7^)} are 
the eigenvectors of A with eigenvalue and {\(3j,Q)} are eigenvectors of B with eigenvalue (3j where (Tj 
and (j shows possible degeneracies. From the assumption of this lemma we know that all of the following 
operators are in J^. 

{\a,,a,){a,,a[\}[j{\P,X,)%X-\} 
Moreover we know that Jif is closed under commutation. So the following operators is also in 

[|a,,cr,)(ai,CTi|, |/3j,Cj)(/3j,Cj|] = ctjI/^j, Cj) cr*) (/3j , | - {f3j,Cj\ai,cr,)\Pj,Q{'^i,'^i\ (9-20) 

From the assumptions of the lemma, we also know that for each pair of a^, Pj there exist some 0"^, such 
that {ai, <Ti\f3j, 7^ and so the following operator is a member of J^. 

e'V., f^.X/S,, 01 - e-^'^l/?,, (9.21) 

where e*^ = {ai,cri\i3j,(j) /\{oii,<Ji\f3j,Cj)\- Also the commutator of this operator with \i3jXj){Pj,Cj \ + 
\(3j , Cj){0j I Cj I is a member of Jif 

= e*«K,a,)(/?,,Cj| + e-^''|/5,,Ci)(a.,a.| + ci|/3,,0)(/?,,Cil + ct|^ 

where ci is a complex number and C2 is real and moreover we have assumed Q =/= Cj. But the terms with 
coefficients ci and C2 are members of Jf and is closed under linear combination. So we deduce that 
e^^\ai,(Ti){(3j,Cj\ + is also a member of Jf. On the other hand, {\f3j,Cj)} is a complete 

basis and so we can expand all {jaijCTi)} in terms of them. This implies that that for all ai,aj and ai,aj, 
the operator \ai,ai){aj,aj\ + \aj,(Tj){ai,cri\ is also a member of J^. So all symmetric operators are in Jif. 
Moreover we know that the commutator of this operator and \ai,ai){ai,ai\ are also members of J^f. This 
would be 

[\ai,ai){aj,a.j\ + \a-j,aj){ai,ai\, \ai,ai){ai,ai\] = \aj,aj){ai,ai\ - \ai,ai){aj,aj\ (9.22) 

So all asymmetric operators are also in J^f and therefore Jif is equivalent to the space of all operators. This 
means that by a sequence of unitary time evolutions commuting with A and unitary evolutions commuting 
with B we can perform all unitaries. 

It was shown in [211 US] that if a compact Lie algebra is generated by {Xi, X2, ■ ■ ■ , X^} then any 
member of the associated Lie group can be generated by a sequence as 6^^*16^*^ ...e'^** where for each of 
these exponentials X is a different member of {Xi,X2, ■ ■ ■ ,Xn}- Moreover the length of this sequence is 
uniformly bounded. Since U (N) is compact we deduce that the length of the sequence we need to generate 
all members of U (N) by a sequence of unitaries commuting with A and unitaries commuting with B is 
uniformly bounded. 
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